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ABSTRACT 

We present imaging and spectroscopic observations from the Interface Region 
Imaging Spectrograph (IRIS) of the evolntion of the flare ribbon in the SOL2014- 
04-18T13:03 M-class flare event, at high spatial resolntion and time cadence. 
These observations reveal small-scale snbstrnctnre within the ribbon, which man¬ 
ifests as coherent qnasi-periodic oscillations in both position and Doppler veloci¬ 
ties. We consider varions alternative explanations for these oscillations, inclnding 
modulation of chromospheric evaporation flows. Among these we find the best 
support for some form of wave localized to the coronal current sheet, such as a 
tearing mode or Kelvin-Helmholtz instability. 

Subject headings: Sun: chromosphere — Sun: flares — Sun: transition region 


Introduction 


The standard picture of the formation of flare ribbons in the sola r chromosphere be¬ 
gins with magnetic reconnection in the corona fiKopp fc Pnenmannlll976l) . As magnetic held 
lines reconnect across a current sheet, coronal hare loops are formed. The plasma along 
these loops is heated to hare temperatures both by the reconnection itself an d by the sub¬ 
sequent contraction of the loop under magnetic tension (ILongcope et al.ll2009l) . The energy 
from this hare pl asma is then transpor ted down the legs of the loop by non-thermal parti¬ 
cles f Brown 1973ll. thermal conduction 1 Craig fc McClvmont 1976 : Forbes et ah 1989 ). wave 
propagation (IRussell fc Fletcheiil2013|l . or some other means, until it reaches the cool, dense 
chromospheric plasma located at the loop footpoints. Energy is rapidly deposited in the chro¬ 
mosphere and transition region (TR), resulting in plasma hows both up and down the loop. 
The uphows are called chromospheric evaporation, and are responsible for hlling the hare loop 
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with hot plasma. Evaporative upflow of hot material is frequently accompanied by down¬ 
ward motions of cooler rnaterial, termed chromospheric condensation flTchimoto fc Kurokawa 
1984J: iBrosius &: Phillips! l2004l: iMilligan fc DennisI 1200911 . Both evaporating and condensing 
components are signihcantly brighter than ambient material, and thus appear as elongated 
emission called flare ribbons. 

As the flare progresses and new magnetic field lines are reconnected, the flare ribbons 
move outward from the po ^rity inversion line, which also marks the approximate location of 
the coronal current sheet flKopp fc Pneuma.nnlll976h . The flare ribbon emission is generally 
considered to represent the footpoints of recently reconnected held lines, and the ribbons are 
therefore an in direct image of the co ronal reconnection process which is imprinted into the 
chromosphere flForbes &: Priestlll984j) . The structure and evolution of the r ibbons can thus 


serve as an observational proxy for ev ents occurring along the current sheet flLongcope et ah 


20071: lQiull2009l: iNishizuka et al.ll2009l) . Signihcant effort has been put into using observations 
of chromospheric ribbons to i nvest i gate the coronal rec onnection (ISchmieder et al.l 119871: 


Falchi et al.l 119971: llsobe et al.N2005 : Mik 


and current sheet electric held flQiu et al 


20021 ). 


enic et al.N2007l) . including the reconnection rate 


The general trend of ribbon spreading during a hare has been established for several 
decades, however it has only been since the early 2000’s that instruments with sufficient 
resolution and cadence, such as the Transition Region and Coronal Explorer ( TRACE, Handy 
et al. 2009), have allowed for observations of small-scale substructure within the ribbon 
itself. Several studies of TRACE observations found that hare ribbons in some events were 
broken into small and tightly-spaced brigh t sources, which were dubbed “compact brigh t 
po ints” (CBPs) ( Warren &: Warshall 2001 : Fletcher fc Warren 2003 : Fletcher et al. 2004). 
In iFletcher et al.l (1200411 . CBPs were found to exhibit a random component to their motion 
in addition to the large-scale ribbon spreading, and they interpreted CBPs as individual 
loop footpoints that wander through the magnetic canopy during the course of the hare. 
Additionally, they found that the brightness of a CBP was correlated to the product of 
footpoint speed and line-of-sight magnetic held, which they interpreted as a measure of the 
local coronal reconnection rate. 

High-resolution observations of coronal loops, meanwhile, have strongly suggested that 
instabilities may occur within the reconnection region which subsequen tly induce oscillations 
in the hare loop flAschwanden et al.lll999l: lOfman fc ThompsonI 1201 111 . One such instability 
which has been studied in the context of the solar corona is the Kelvi n-Helmholtz instability , 
which occurs at a huid interface with a discontinuity in how speeds (lUchimoto et al.lll99111 . 
KH in stabilities have been specih cally invoked to explain observed oscillati ons in coronal 
loops flOfman &: ThompsonI 120111) and auroral spirals in the magnetopause flLvsak fc Song 
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1996|) . Another instability that has attracted interest is the tearing mode (TM) instability, in 
which magnetic isl ands spontaneously form and grow in a current sheet during reconnection 
flFurth et ahlll963l ). The TM instability is believed to be important in permitting the fast 
reconnection rates r equired for s olar flare energy release in so-called “impulsive bursty re¬ 
connection” models fjPriestlll985l) . and numerical simulations of reconnection h ave frequently 


invoked some forrn of ve locity shear or TM instability along the current sheet flKarpen et ah 
I 995 I: iKliem et ahlbnooh . 


Despite some successes flAschwanden et ahl Il999l: lOfman fc ThompsonI I 2 OIII) , it has 
proven difficult to unambiguously resolve either current sheet instabilities or the associ¬ 
ated flare loop oscillations, due in part to lower emission intensities and line-of-sight effects 
in the corona. The chromospheric flare ribbons, on the other hand, are much more intense 
relative to the background emission, and any footpoint brightenings or ribbon motions are 
more easily distinguished from noise. The connectivity between the reconnection region and 
the chromospheric footpoints means that a global instability in the coronal current sheet 
might be expected to imprint itself into the ribbon evolution. Further, since the entire cur¬ 
rent sheet is presumably connected to the ribbons via the reconnecting flare loops, a current 
sheet instability would likely manifest as a substructure in the ribbons which could be easily 
distinguished from random CBP motion by a coherent pattern. This assumes, of course, 
that held-aligned transport mechanisms translate energy consistently down the loops. An 
alternative is that the observed ribbon patterns might be a result of the transport rather 
than the source. The conventional understanding, based on numerous observations of coher¬ 
ent ribbon motion at two magnetically conjugate sites, is that the transport does not play a 
signihcant role. 

In fact, the Interface Region Imaging Spectrograph (IRIS, De Pontieu et ah 2014), which 
makes high spatial and temporal resolution observations of the chromosphere and TR, has 
captured instances of such coherent sub s tructu re in flare ribbons. One example is a recent 
study of an X-class flare bv lhi fc Zhand (120151) . where the authors identify a quasi-periodic 
slipping motion of flare loop footpoints which is observed as a series of small bright knots in 
Si IV which oscillate along the ribbon. They interpret this slipping behavior as an apparent 
motion of footpoints, which brighten and fade as quasi-periodic slipping reconnection drives 
impulsive footpoint heating. Finally, based on the timing of the oscillations, they suggest 
that the slipping reconnection may be due to varying densities in the current sheet which 
are driven by p-mode oscillations above the sunspots. 

In this paper, we present an analysis of an IRIS observation from 2014 April 18 of a 
two-ribbon flare which displays coherent substructu re of the ribbon d uring the impulsive 
phase of the flare, very similar to that reported by Li &: Zhanel ( 2015 1. We hnd support 
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in the data for a scenario of a flare loop undergoing elliptical oscillations, driven by either 
a TM or KH instability in the current sheet. Our paper is outlined as follows: in Section 
m we describe the data and the details of the flare event. Then, in Section [31 we describe 
the evolution of the flare ribbons and our method for spectral line htting, determine the 
Doppler velocities for Si iv within the ribbon, and lastly discuss both the behavior of other 
spectral lines and compare the substructure in both ribbons between conjugate points. Next, 
in Section 01 we consider and dismiss two alternative scenarios for generating the observed 
ribbon substructure before proposing our scenario for instability-driven elliptical oscillation. 
Finally, in Section |5l we discuss some of the possible implications for our results. 


2. Observation 


2.1. Instrument 


The Interface Region Imaging Spectrograph (IRIS) is a space-based observatory in low- 
Earth orbit that was launched on 2013 June 27. The primary instrument onboard IRIS 
is a dual-range UV spectrograph (SG) with 0".16 pixels and an effective effective spatial 
resolution of 0".4. The SG slit is 0".33 wide and 175" long, and covers FUV passbands 
from 1332 A to 1358 A and 1389 A to 1407 A and an NUV passband from 2783 A to 
2835 A. These passbands include lines formed over a wide range of temperatures from the 
photosphere (5000 K) to the corona (1 to 10 million K). Two hducials (dark bands where 
light is excluded from the SG) are located along the slit which provide alignment and spatial 
conte xt, and the entire SG c an be rastered back-and-forth to provide coverage over a 2-D 
area (jPe Pontieu et al.ll2014j ). 


The other instrument featured onboard IRIS is a slit-jaw imager (SJI) which records 
context images of the observation region on either side of the SG. The IRIS SJI includes four 
wavelength passbands, including two transition region lines (Si iv 1400 A and G ll 1335 A), 
with a £eld-of-view (FOV) of 175" x 175", and t ypically operates at one -third the cadence of 
the SG (due to cycling through the passbands) fjPe Pontieu et al.ll2014l) . The inclusion of a 
SJI, with the same FOV, resolution, and alignment as the SG allows for detailed contextual 
knowledge of the slit placement relative to observational features such as sunspots, ribbons, 
and coronal loops. It also facilitates very precise co-alignment (< 1") between IRIS and 
other instruments, such as SDO/AIA, by aligning visual features in the SJI to images from 
those instruments. 
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2.2. Flare 

The event investigated in this paper is a GOES M7.3-class flare that occnrred on 2014 
April 18. The flare was located in NOAA Active Region 12036, which on that day was 
located at approximately 500" West and 200" Sonth from solar disk center. The flare start 
time was 12:31 UT and the GOES X-ray flux peaked at approximately 13:03 UT, with the 
event mostly concluded by 13:20 UT. The large-scale structure consists of two chromospheric 
ribbons, one on the east side and one on the west side of the flare, visible as strands of intense 
1600 A emission at 12:46:40 UT in Figure [1] (note that the color scale is a reversed log^g 
black-white (RLBW), with white as lowest intensity and black as highest intensity). We 
have labeled the east and west ribbons accordingly as “ER” and “WR”, and we refer to the 
ribbons as such throughout the text. Both ribbons generally follow a single distinct path 
from southeast to northwest, although the ER is more strongly tilted toward a north-south 
orientation. There are also occasional offshoots and patchy regions of intensity which evolve 
with time. 

One notable exception to the general SE-to-NW single-track ribbon paths can be seen in 
the WR starting at around 550" West and 200" South (Figured]). At this point, moving SE 
to NW, the ribbon clearly splits into two channels which merge back together farther to the 
NW, resulting in an isolated island between the ribbon channels. This split ribbon feature 
persists over most of the evolution of the flare and appears in both SDO/AIA 1600 A and 
IRIS SJl 1400 A, but does not appear in the SDO/AIA 171 A passband. The most likely 
origin of this feature can be seen in the SDO/HMl magnetogram for this region taken in the 
middle of the flare duration (12:58:19 UT) as shown in Figure [2l where positive and negative 
LOS polarity are colored white and black, respectively, and the ER and WR positions are 
traced for context as the white and black lines. In the center of the split WR island we 
observe a significant neutral gap running diagonally through the positive magnetic region 
centered at 575" West and 200" South. The split 1600 A ribbon appears to track around 
either side of this neutral channel. 

There is also a prominent arcade of coronal loops which appears in the SDO/AIA 171 A 
passband approximately 30 minutes after the first appearance of the flare ribbons. Figure [3] 
(also RLBW) shows a 171 A image of the flare region taken at 13:12:37 UT showing the flare 
ribbons shortly after they first appear. For comparison with Figure dl we have traced the 
position of the ribbons as jagged white lines. The loops visible in Figure [3] all follow roughly 
E to W paths between the two ribbons, although the exact loop connections between the 
ER and WR are not easy to establish because of the time delay for the appearance of the 
loops in this passband. We also searched for loops at earlier times using higher temperature 
SDO/AIA passbands, such as 131 A and 193 A, but saturation and lack of sharpness in 
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these passbands relative to 171 A impedes clear identification. Our main use for the AIA 
imaging data is to identify magnetic conjugacy. We use 171 A images for this because they 
are sharpest and will probably extend to the lowest points along the loops. The 171 A image 
clearly indicates a dipolar component to the magnetic geometry, which we confirm with the 
SDO/HMI magnetogram in Figure [2l We clearly see two dominant polarities, negative to the 
east and positive to the west, which the ribbons track quite closely and which are connected 
by the flare loops. Finally, we note there was a coronal mass ejection (CME) associated with 
this flare (observ ed by LASCO'); the flux rope eruption responsible for the CME is discussed 
in more detail in ICheng et al.l (120151) . 


The IRIS SG data we use for this study is drawn from a sit-and-stare observation (i.e. 
no rastering of the SG slit) which began at 12:33:38 UT, shortly after the GOES start time 
for the flare, and which continued until 17:18:10 UT. This observation thus covers nearly 
all of the impulsive rise phase and the entire decay phase of the flare. The SG cadence for 
this observation was ~9 s, while the SJI A cadence for 1400 A was ~27 s, and the SJI FOV 
was initially centered at 550" West, 230" South (note that the spatial metric at this position 
on the disk is ~900 km/arcsec). An example SJI 1400 A image, taken at 12:46:34 UT 
(approximately the same time as Figured]), is shown in Figured] (we again use RLBW). For 
context, note that the IRIS SJI FOV is shown as a black inset box in Figures [IlEl and that 
the IRIS SG slit position is indicated in those three figures by a vertical dashed black line. 
Finally, note that the IRIS SG slit is conveniently positioned across the WR, directly over 
the location of the magnetic gap and ribbon island mentioned above. In summary, the slit 
placement, high cadence, and consistently stable rotation tracking during this observation 
make it ideal for a spectroscopic study of the fine structure of flare ribbons and its evolution. 


2.3. Wavelength correction 


The rest wavelengths of spectral lines for the IRIS SG are generally dynamically shifted 
to some degree due to the orbital motion and thermal variations of the IRIS spacecraft and 
must be corrected for every spectrograph exposure. To do this we use the SSWIDL routine 
iris_orbitvar_corr_12, which fits the Ni I 2799.474 A line to determine the appropriate 
wavelength shifts for the FUV bandpass (iDe Pontieu et al.ll2014i) . Even after applying this 
correction, however, we find that the peak of the O l 1355.598 A is systematically redshifted 
during the observation by 01 A, even though as a photospheric line it should generally 
be stationary. We therefore subtract this additional shift from the FUV wavelength axis to 
make the O l line sta tionary, and subsequentl y reference all other lines to their GHIANTI 
database line centers fjiRIS Technical Note 20l) . 
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Results 


3.1. Ribbon evolution 


The flare ribbon appears as an intense band in the northern portion of the IRIS SJI 
beginning around 12:46:34 UT. The overall orientation of the WR is from SE to NW, and in 
the vicinity of the SG slit the ribbon is running very nearly due E-W, across and perpendic¬ 
ular to the slit. As the WR evolves it slowly drifts S, generally maintaining its perpendicular 
orientation to the slit. This motion is away from the polarity inversion line (PIT), in agree¬ 
ment with the classic picture of ribbon spreading (iKopp fc PnenmaunlllQTbh . Figure [5] shows 
a series of 16 images from the SJI 1400 A passband, taken from the black inset box near 
the top of Figure m at ~1 minutes intervals starting at 12:46:34 UT (the same time as Fig¬ 
ure S]) and continuing to 13:00:41 UT. The color scale is RLBW. The inset frame has been 
chosen to cover the area of the SJI where the northern branch of the WR crosses the SG 
slit. The reader should note that the time axis for these 16 frames runs from left-to-right 
in the top row, then right-to-left in the second row, and so on row-to-row as indicated by 
the arrows to form a movie of the ribbon evolution. The SG slit can be seen as a pale line 
running down the middle of each frame, with the upper fiducial just visible at the bottom of 
each frame. Intermittent saturated pixels appear to the left of the slit in some frames (for 
example, 12:52:41 UT). 

In addition to the large-scale evolution of the ribbon shown in Figure |5l we also note 
a distinct substructure that evolves on a smaller scale and at a faster rate than the overall 
ribbon motion. This substructure appears in most frames as a jagged sawtooth pattern that 
cuts across the slit, with multiple patches of bright emission orient ed diagonally to th e slit 
that break up the ribbon. A very similar structure was observed by iLi &: Zhand fl2015h and 


called “slipping reconnection” by them. Since this term also refers to a model flAulanier et ah 


2006|), which may or may not be pertinent to the observation, we use the term sawtooth here 
instead. Times at which this pattern appears most prominent include 12:50:49 UT, 12:52:14 
UT, and 12:55:32 UT. An inspection of this sawtooth pattern from frame-to-frame reveals 
that the diagonal features appear to slide across the slit from east to west as the ribbon 
drifts south over time. Our best estimate of its pattern speed is Vgt ~ 15 km s“^, parallel to 
the ribbon itself. One example where this is most readily apparent occurs between 12:53:38 
and 12:55:32 UT. In these three frames a diagonal sawtooth begins just to the left of the 
slit, moves west such that one minute later it is crossing the slit, and one minute after that 
is seen on the right side of the slit. 

Since the sawtooth pattern in the SJI 1400 A is moving across the SG slit, we would 
expect the pattern to also be reflected in spectral intensity for the Si iv 1403 A passband. 
















To calculate the total intensity in this passband, we hrst divide the data for each SG frame 
by the appropriate exposure time to obtain normalized units of DN s“^. We then reset all 
negative values in the passband (which occur near the passband edges) to 0, and finally sum 
the data over the entire Si IV 1403 A passband. In the upper plot of Figure [HI we show 
the resulting time-distance stackplot for the time range 12:45 UT to 13:05 UT. The y-axis 
shows the heliographic position of the slit pixels at at the beginning of the observation; 
there is some wobble to ^/-pointing of the instrument, but it is less than 0.3" over the entire 
observation. Data values vary between 20 DN s“^ and 64000 DN s“^. We again note the 
southward motion of the ribbon, beginning at ~184" S and subsequently moving south to 
~188" S. We also observe immediately that the sawtooth pattern is indeed present in the 
SG data starting around 12:48 UT and ending at around 13:02 UT, and that three teeth 
in particular can be clearly seen between 12:51 UT and 13:01 UT. There is also a feature 
that moves up from the south beginning at around 12:47 UT, however this is a piece of the 
northward-moving southern branch and not part of the ribbon sawtooth structure. 

In the lower plot of Figure [6l we have replotted the intensity stackplot (with the time 
axis in seconds this time), and we have overlaid a red contour line around the section of the 
ribbon emission for which the sawtooth pattern is the clearest. Henceforth, when referring 
to the “sawtooth”, we will be referring to this outlined region of pixels. The numbers 
1-6 indicate six peaks in the sawtooth oscillation, with dashed lines also indicating the 
approximate beginning time of a southward motion in the sawtooth. The oscillation period 
varies between 80 and 190 s, averaging ~140 s, consistent with sawtooth structures ~2 Mm 
long moving across the slit at ~15 km s“^. We have also overlaid the sawtooth centroid 
position (as a function of time) as the blue line, as well as a linear best fit to the sawtooth 
mean as the orange line. The linear fit has a southward velocity of 1.6 km/sec, while the mean 
position moves at velocities ranging from ±20 km/s. This confirms the much faster motion 
of the sawtooth substructure relative to the overall ribbon motion. Finally, we observe that 
the shallow sides of the sawtooth pattern are all directed away from the PIL located between 
the ER and the WR. 


3.2. Spectral line fitting 


The WR sawtooth substructure is seen in the 1400 A SJI passband. These images 
are dominated by the two Si iv lines (139 4 A and 1403 A) which are recorded in the IRIS 
SG FUV 2 range fjPe Pontieu et al.ll2014l) . These two lines are ideal for examining flare 
ribbon plasma, since at ~80,000 K they are located within the transition region at the 
footpoints of the flaring loops. Additionally, they are typically strong lines that dominate 
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their respective spectral regions and therefore are easily distinguished from the lines of other 
ion species. In Figure [7] we have plotted the exposure-normalized data for the two lines for 
several representative positions and times in and around the sawtooth pattern: 1394 A as 
asterisks and 1403 A as squares. The x-axes are in units of km s“^ with n = 0 referring to 
the rest wavelength defined as described above, where redshifts and blueshifts are defined as 
positive and negative velocities respectively. The y-axis is arbitrarily scaled for each pixel, 
and the time appears to the left of each plot. At the upper right is the outline of the sawtooth 
from Figure [HI with symbols arranged in five rows “A”-“E”, corresponding to a row of 
spectral plots, which indicate the position and time of each spectral plot. Finally, note that 
the 1394 A data have all been scaled by a factor of 0.5, so that the two spectral lines will 
appear on the same scale (and see the hnal paragraph of this section for more detail on this 
choice). 


The positions and times (henceforth “pixels”) chosen for Figure [7] have been chosen to 
illustrate several general features of the Si iv spectral lines in this event. First, the majority 
of the pixels within the sawtooth appear to contain two Gaussian components, as first noted 
by ICheng et ahl (120151) in the early ribbon development. The two components are present 
for both the 1394 A and the 1403 A lines, clearly demonstrating that both components are 
from Si IV and not from an accidental blend with another spectral line. A fraction of pixels 
late in the sawtooth contain one or more additional components; however, these occur only 
in a minority of cases and typically have small amplitude compared to the two dominant 
components. The lack of additional components, or of significant non-Gaussian tails in 
the spectral lines, suggests that the IRIS SG is not observing many separate loop footpoints 
within a single pixel, but rather only one or two distinct footpoints. We also observe that each 
component appears to persist and to evolve within its row. Both components are Doppler- 
shifted, with one component consistently redshifted and the other switching between redshift 
and blueshift, and both the Doppler shift and the velocity separation between the components 
evolves with time. Finally, we note that neither component consistently dominates, and for 
some pixels they have comparable magnitude. 


Within the sawtooth, there are numerous cases where at least one, and sometimes both, 
of the Si IV lines saturate the IRIS SG pixels. This is true for ~17% of pixels for 1394 
A and ~5% of pixels for 1403 A, and makes accurately fitting a two-component Gaussian 
problematic. Except for saturation, however, most positions display behavior similar to that 
discussed above, although we discuss some additional discrepancies below. We therefore 
assume that a two-component Gaussian will be appropriate for the majority of the pixels 
within the sawtooth. 


Since there are too many positions and times within the sawtooth to manually supervise 
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the fitting, we have developed an automated procedure for eliminating poor fits. We note 
hrst that the chi-squared measure for the hts were found not to neatly discriminate between 
fits which were deemed good and bad by visual inspection. Therefore, in our procedure the 
two-component £t is rejected if any one of the following criteria are met for either component: 


1. any pixel in the band is greater than 2000 DN s“^ (approximate saturation limit); 

2. the Gaussian amplitude is less than 10 DN s“^; 

3. the line center is greater than ±0.625 A (~134 km s“^) from the rest wavelength; 

4. the Gaussian width is greater than 0.625 A. 


We find that ~21% of 1394 A and ~7% of 1403 A attempted fits are rejected using this 
method; roughly two-thirds of these rejections are due to condition 1, while the rest are 
from the other three conditions. Because of the greater percentage of successful hts for 1403 
A, we will use only that spectral line for the remainder of the analysis starting in Section 
13.31 In a small number (<1%) of cases, only one component is the cause of the rejected fit; 
however, due to the rarity of these cases, we simply reject both components. Finally, the two 
components are sorted according to their Doppler shifts, and we label the bluer component 
Cb and the redder component Cr. 


In Figure[7]we have plotted the resulting two-component Gaussian hts for the 28 selected 
pixels over the data for those pixels, for both the 1394 A (as the solid orange curve) and 
the 1403 A (as the solid green curve) spectral lines. The dashed lines correspond to the 
individual components Cb and Cr for each of the Si IV spectral lines, with cyan/violet used 
for the respective 1394 A spectral components and blue/red used for the respective 1403 A 
spectral components. We note that our criteria are generally adequate to reject bad hts, and 
the hts that remain conform to the data for both spectral lines quite well. We also note that 
the ht parameters of the two components are generally similar for both spectral lines, which 
adds conhdence that the htting procedure is hnding the same two major components. We 
hnd that Cr is generally narrower than Cr (for ~88% of pixels). The average non-thermal 
widths, given by = \/o'^ ~ o'tb ~ where crth = 6.86 km s“^ and dinst = 3.9 km s“^ 
(jPe Pontieu et al.N2015l) . for the two components are a® = 19 km s“^ and cr^ = 38 km 




s^". We note that the non-thermal widths for Cb are quite consistent with those found in 
quiescent active regions (jPe Pontieu et al.ll2ni5h whereas Cr shows much greater broadening, 
but we will defer discussion of this fact until later. 


In an optically-thin plasma w ithout geometrical effects , the 1394 A and 1403 A Si iv 
lines will have a 2:1 intensity ratio (IMathioudakis et al.lll999h . The majority of the sawtooth 
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pixels display this 2:1 ratio: e.g. the majority of the scaled data for rows A, D, and E all lie 
on top of each other. Some pixels, however, depart from the strict 2:1 line ratio, snggesting 
some optical depth effects. This departnre is more common in Cb, especially when it is very 
narrow, as in rows B and C. It can also occnr in Cb, as it does in row D. In some cases the line 
ratio is less than two, while in other it is greater (as in Cr of row D). The former sitnation 
is co mmonly attribnted to loss fr om the thicker line (1394 A) by scattering ont of the line of 
sight flMathiondakis et al.l 1199911 . The latter c an occnr in more complicated geometry when 
photons are scattered into the line of sight ( Kerr et ah 20051) . The fact that we observe 
both cases in neighboring pixels snggests, not snrprisingly, that the emitting plasma has a 
complex geometrical strnctnre over scales comparable to the IRIS resointion. Almost all 
optical-depths effects we observe are similar to those shown in rows B, C and D. The line 
ratio remains very close to the optically thin valne of 2:1, and the thicker component (1394 
A) is still well £t by a Ganssian. We therefore conclnde that Si iv is never far from being 
optically thin and that Doppler shifts and line width are well measnred by the Ganssian fit. 


3.3. Doppler velocities 

In FignrelHl we plot the Doppler velocities determined above for both components nsing 
a bine-red color scale on time-distance stackplots {Cb in the npper panel and Cr in the lower 
panel). The time axes follow the same labeling scheme as in Figure El and the blue-red color 
range runs between ±40 km s“^ for Cr and between ±80 km s“^ for Cr (in this system, we 
define redshifts as positive and blueshifts as negative). We have also outlined the sawtooth 
pattern from Figure El with the solid black line to help guide the reader in determining the 
relative position of features. Finally, we have color-coded any pixels not £t (for any reason) 
by the above automatic routine as gray. 

There are several important features to note in the Doppler velocity data, mostly for 
component Cr- The first is that, as mentioned above, the Doppler shifts for component 
Cr are sometimes red and sometimes blue. However, we now observe that these shifts are 
not random, but are clustered into large regions of redshift and blueshift. Additionally, we 
note that the blueshifts tend to be concentrated near the upswings and peaks in the saw¬ 
tooth pattern, whereas the redshifts tend to be on the downswings and troughs (particularly 
prominent in the teeth marked “3” and “4”). Although it is not as apparent, this same 
pattern is reflected in component Cr, except that the redshift is merely weakened instead of 
shifted to the blue. 

A different method of visualizing the velocity evolution is to plot the Doppler velocities 
for all positions in a single time slice of the sawtooth as a function of time, as we have 
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in Figure |9l Components Cb and Cr are displayed as the black squares and asterisks, 
respectively, and the mean velocity within each column is plotted as the solid blue (for Cr) 
and red (for Cr) lines. In this format we can clearly observe that, starting at ~1100 s, there 
begins a distinct red-blue oscillation in (7^, apparent in the mean velocity, which continues 
for at least three full periods of 100-200 s and amplitude ± 20 km s“^. Finally, we note 
that there is a similar oscillation in the redshifted Cr Doppler velocity, which is noisier than 
for Cb but still clearly at the same frequency as and in phase with the Cr oscillation. The 
similar frequency and phase of these Doppler oscillations strongly suggests that both Si iv 
components are being shifted by a common motion of the plasma. 


3.4. Additional spectral lines 


There are two additional pairs of spectral lines included in the IRIS SG that are of use 
for observing the upper chromosphere and TR. One pair is the O iv 1400 A and 1401 A lines 
included in the FUV 2 passband. These lines are formed higher in the TR than the Si iv 
lines, at ~150,000 K, and as f orbidden lines they are often of use for determining the plasma 
density at that temperature flFlower fc NussbaumerlllOTSl) . We hnd that the time-distance 
stackplot of the O iv 1401 A total line intensity is virtually indistinguishable in structure 
from that of Si IV 1403 A, with the sawtooth pattern appearing at an identical position and 
time. Unlike the Si iv prohles discussed above, however, neither of the O iv lines have any 
obvious two-component structure and appear generally well-ht by a single Gaussian. This 
may be because the O iv lines have count rates signihcantly lower than the Si iv lines, or 
it may be due to the specihcs of the formation mechanism for the forbidden O IV lines. 
Although we do not discuss it further in this work, we hnd that the single component O IV 
1401 A Doppler shifts generally follow that of component Cr of Si iv 1403 A, with similar 
magnitudes of blueshifts and redshifts located in the same parts of the sawtooth. 

The other pair of lines are the G ll 1335 A and 1336 A lines in the FUV 1 passband. 
These lines are formed in the upper chromosphere, at ~20,000 K, and are thus expected to 
be (somewhat) closer to the footpoint than Si iv and generally will not display the same 
dynamics. However, we hnd that G ll 1335 A and 1336 A lines also display nearly identical 
behavior to Si iv, with the sawtooth intensity pattern appearing over the same spacial and 
temporal scales. We also hnd t hat b oth G ll lines possess two-component Gaussian prohles, 
as was noted in iGheng et al.l (120151 ). and although we did not perform the same detailed 
analysis of the Doppler shifts we note that these components display very similar shifts to 
the Si IV lines at a several selected pixels within the sawtooth. 








13 


3.5. The East ribbon 


Although the IRIS SG FOV only allows for spectral analysis of the WR, we note that 
the northern end of the ER was located with the IRIS SJI 1400 A FOV (Figure 0]). From the 
171 A SDO/AIA image in Figure [3l we observe a set of flare loops that emerge from the WR 
at the IRIS SG slit location which curve first North and then South and West to end at the 
ER at the location indicated in that figure by a vertical black “I”-shaped line crossing the 
ER. We extract an artificial SJI slit, one pixel wide and located along the vertical black “I”- 
shaped lines in Figures [T] & 01 to investigate the ER intensity behavior in the same manner 
as the WR. The resulting intensity time-distance stackplot is shown in the middle panel of 
Figure [TUI with the original SG stackplot from Figure 0] plotted in the bottom panel. 

We note immediately that there are two diagonal bands that begin at very nearly the 
same times as the peaks marked “3” and “4” in the original sawtooth. Of course, the lower 
time resolution of the SJI (~27 s) means that any features in time are somewhat less clear 
than for the ~9 s cadence SG. However, the SJI bands both begin within 60 s of their 
companion SG sawtooth features, and their duration (120-180 s) and spatial extent (2-3") 
are nearly identical. We also note that the slow motion of the band is also directed away 
from the PIL, as it was for the SG sawtooth, indicating that the phase of these features is 
the same. 


We repeat this procedure using the co- aligned AIA 1600 d ata, which is often dominated 
by the G IV lines at 1548 A and 1550 A f Lemen et al.l 12012 ). These lines are formed at 
~100,000 K, slightly hotter than for Si iv. The resulting intensity stackplot is shown in the 
top panel of Figure (TUI Once again, we note that the lower spatial (0".5) and temporal (24 
s) resolutions smear out the features, but we can still make out the same diagonal bands 
marked “3” and “4” at approximately the same positions and times as before. Obviously, we 
cannot be absolutely certain that these ER footpoints are conjugate to the WR footpoints 
in the IRIS SG. We also note that the bands present in the ER are not completely identical 
to the sawtooth in the WT, and that the oscillations “1”, “2”, “5”, and “6” from the WR 
sawtooth are conspicuously missing in the ER. Further, those oscillations do not appear for 
any other artificial slit positions on either side of the position shown, and as the artificial slit 
is moved away from that position the two oscillations that do appear are not as apparent. 
This indicates some degree of asymmetry between the two ribbons. However, we are fairly 
confidant stating that the conjugate behavior of the intensity stackplots indicates that both 
ribbons are experiencing the same oscillation, and that the close relationship we observe 
between the ER and WR strongly implicates a source somewhere high in the corona. W e 
therefore rule out any local mechanism, such as p-mode oscillations (iLi fc Zhang! l2015h . 
which would not be linked by the coronal magnetic field. 
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4. Interpretation 

In the upper panel of Figure fTTl we have plotted the mean Cb Doppler velocity in km 
s“^, determined in Section [3731 and shown in Figure [9l against the relative difference in arcsec 
of the mean sawtooth position from the linear £t (effectively the orange line subtracted from 
the blue line in Figure |H]). The result is a phase portrait for the sawtooth ribbon oscillations. 
Note that the actual sawtooth data are represented in the plot by the colored squares, with 
the connecting colored lines added to guide the reader. Time is encoded by the color of 
the squares and connecting lines starting with purple at the beginning of the sawtooth and 
ending with red, as shown by the color bar on the right side of the figure. 

The central feature of this phase portrait is that, despite some meandering of the data, 
many of the data appear to follow a diagonal trend that runs from upper left to lower right. 
The solid black line represents a linear fit to the full data set within the sawtooth, and 
gives an overall slope of ~14 s (using ~900 km/arcsec). However, we note that the diagonal 
trend becomes more apparent when we consider only individual pieces of the sawtooth. In 
particular, in the lower panel of Figure [TTl we have grayed-out all the data except that for 
sawtooth “4” in Figures [6] and [8] (beginning at ~1400 s and ending at ~1520 s). These 
data are clearly aligned along a diagonal trend, and the dash-dotted line represents a linear 
fit to just these data with a slope of ~20 s. In the context of the discussion so far, the 
diagonal trend makes sense: the strongest blue and red Doppler shifts tend to be concentrated 
respectively at the peaks and troughs of the sawtooth pattern. The phase portrait presented 
in Figure ITT] therefore reveals a critical property of the ribbon substructure. Namely, that 
the LOS Doppler velocity and the position of the sawtooth pattern are 180° out-of-phase. 

Before presenting our proposed scenario for the generation of the sawtooth pattern, we 
first consider two alternative scenarios that would be probably serve as the most obvious 
first candidates. The first scenario is a simple harmonic oscillation (SHO), where the loop 
is waving forward and back within a plane. This could come from a linearly polarized loop 
oscillation and would create redshifts and blueshifts in the loop plasma as the loop motion 
carries it toward or away from the observer. However, these Doppler shifts would occur only 
as the loop was moving from North to South or from South to North, and cannot account 
for the Doppler shifts at the peaks and troughs of the wave. We also note that the phase 
portrait for an SHO is an ellipse in position-velocity space, resulting from the 90° phase 
difference between position and velocity. The SHO scenario is therefore inconsistent with 
the results from the ribbon sawtooth. 

The second scenario which we consider is also an SHO ribbon oscillation, but instead of 
the loop motion creating the Doppler shifts it simply modulates pre-existing plasma velocities 
within the loop. In this scenario, the loop tilts first toward and then away from the observer. 
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the Doppler shifts of plasma flowing along the loop would be changed from red to blue 
(or vice versa, depending on the flow direction). Assuming the correct conditions, this 
modulation could provide the correct 180° phase difference. It is readily seen, however, that 
this scenario requires the loop oscillation to be at least partly on the far side of 90° inclination 
from the LOS, in order to switch the sign of the Doppler shift. From the 45-minute HMI 
vector magnetogram at 12:58:19 UT we find that the magnetic field in the sawtooth region 
is angled only ~40° from the LOS. This would require a substantial (>50°) oscillation in 
order to switch the Doppler shifts, which does not seem reasonable given the observations. 

The 180° phase difference between the LOS velocity and the apparent wave position 
suggests that the flare loops (and the fluid elements of plasma along them) are undergoing 
motion that is alternately parallel and perpendicular to the LOS. One example of a class of 
waves which produces such motion are surface waves, such as occurs in deep water driven 
by wind. In a surface wave, traveling in a direction x perpendicular to the surface normal z 
with no mean flow, the displacement of a fluid element hx is given by 

(x, z, t) = [— sin {kx — ut) + cos {kx — ut) e ^]; (1) 


and the the velocity v is given by 

V (x, z, t) = [cos {kx — ojt) bx -|- sin {kx — ujt) e^], 


( 2 ) 


where A is the wav e amplitude, u is the wave frequency, k is the wave number, and ej is 
the unit vector fjPhillipsI 119771) . As can be seen from these equations, a fluid element in 
this wave will trace out an ellipse in (x, z) as the wave passes, with a 90° phase difference 
between the two velocity components and a total 180° phase difference between a given 
velocity component and the position component perpendicular to it. 

We also recall that the symmetry between the two ribbons suggests a coronal source to 
the wave. It therefore seems likely that an instability in the current sheet during reconnection 
is driving the sawtooth. Two such instabilities which have been shown to occur in coronal 
current sheets during reconnection, and which ca n lead to oscillations in the subsequent flare 
loops , are the Kelvin-Helm holtz (KH) instability fjUchimoto et al.lll99ll:l0fman fc Thompson 


2011 : Foullon et ah 2011 ) and the tearing-mode (TM) instability f Furth et ah 19631) 


Like 

the water wave described above, both of these plasma instabilities are confined to a surface 
(the current sheet) and produce elliptical motions. In the upper panel of Figure [121 we 
present a schematic conception of the scenario by which these instabilities could produce the 
observed features. Looking north along the flare arcade, the LOS comes roughly from the 
upper left. The outer two black lines represent a pair of reconnecting held lines, where the 
reconnection is occurring at the star along the current sheet (vertical dashed line). We then 
hypothesize that either a KH or TM instability occurs, which results in elliptical motion of 
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the field lines during their reconnection (indicated by the arrowed circles). The oscillation 
frequency is low enough that the loop responds rigidly, as a whole, pivoting at some depth 
below the chromosphere. Finally, the inner black line shows a previously reconnected loop 
that has contracted and then cooled to become visible in 171 A. By this point the loop is 
disconnected from the current sheet and oscillation has ceased. The instability drives the 
elliptical oscillation, but also appears to have a phase velocity within the sheet, producing 
the pattern motion at Vgt — 15 km s“^. 

Meanwhile, at the loop footpoint we see the effects of the instability-driven oscillations 
on the ribbon. The dotted box in the upper panel of Figure [12] indicates the zoom-in region, 
which is shown in the lower panel of Figure [T2j Here, the oscillating loop is shown in two 
positions as the solid arrowed lines, with the LOS again from upper left as indicated. The loop 
is shown anchored somewhere deeper in the atmosphere, and the elliptical loop oscillation is 
indicated by the ellipse at the center and by the ® and © respectively representing motion 
into and out of the page. The observed amplitude of the oscillation will be the diameter 
of the ellipse; from above we find this to be ~l-3" or ~0.7-2.2 Mm. We note that the 
only motion of the loop that results in Doppler shifts is that of the loop moving along the 
LOS, as indicated by the blue and red arrows, and that the loop motion at the © and 0 is 
perpendicular to the LOS and generates no Doppler shifts. Therefore, as we observe, the blue 
and red Doppler shifts will occur at the northern and southern extremes of the oscillation, as 
in Figure [HI and that the phase difference between the apparent position and Doppler shift 
will be 180°. 


Obviously, in order to produce the observed Doppler shifts in Cb of ±20 km s“\ the 
loop must be rotating around the circle indicated in Figure [12] at approximately that same 
speed (technically it would need to be somewhat faster, due to the inclined LOS). Given 
that the observed amplitude of the oscillation (equivalent to the diameter of the circle) is 
~0.7-2.2 Mm, we note that a rotation speed of ~30 km s“^ (allowing for LOS effects) would 
give an oscillation period of 70-230 s. This is very similar to the observed range for the 
period of the sawtooth oscillation, which varied between 80-190 s. 


Finally, we believe that the two components Cr and Cb occur in distinct sets of loop 
footpoints, which are too closely spaced to be resolved by IRIS. Both sets of loops are partici¬ 
pating in the elliptical oscillation generated by the instability described above. However, the 
greater non-thermal broadening of Cr, as well as the nonzero average redshifts, suggests that 
those footpoints are undergoing chromospheric condensation, presumably as that plasma is 
being directly energized by reconnection in the instability region. The Cr component, on the 
other hand, has non-the rmal broadening consistent with previous quiescent AR observations 
(jPe Pontieu et al.[[2015fl . and displays an approximately zero average Doppler shift during 
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the oscillation. This suggests that this set of footpoints have not been energized into conden¬ 
sation by coronal reconnection, even though those loops are still being elliptically oscillated 
by the instability. We recall, however, that the Cb component is still signihcantly bright¬ 
ened during the oscillation; it is possible that this plasma is being energize d into enhanced 
emission via energy transfer from the loop waves flRussell fc Fletcher! l2013l) . 


5. Discussion 


In this paper, we have presented an analysis of a two-ribbon flare using IRIS and AIA 
imaging and spectral observations. We found, along with the usual ribbon spreading at 
~l-2 km s“^, that the ribbons show a distinct sawtooth substructure with a scale of 1-2 
Mm. The sawtooth appears in both ribbons, maintains its shape coherently over time, and 
drifts along the ribbon from east to west at a speed of Vst — 15 km s“^. The period of 
this sawtooth oscillation was found to average ~140 s, and the oscillation speed is ~20 km 
s“^. We also found that the sawtooth substructure appears in a variety of spectral lines, 
including Si IV, O iv, and C ll, with nearly identical amplitude and phase across a wide 
range of temperatures. We observed that the line prohles for Si IV and C l l in an d around 
the sawtooth have two major components (previously noted bv ICheng et al.l (120151) h and we 
developed an automated routine for fitting the two Si iv lines observed by the IRIS SG. O 
IV does not clearly show two components, but we noted that this may be due to the weaker 
intensity or different formation mechanism of the forbidden O IV lines. We showed that the 
two components of Si iv tend to persist over time and we identified the two components by 
their relative Doppler shifts, with the redder component labeled Cr and the bluer component 
Cb. 


We found that Cr was redshifted by an average ~50 km s“^ during the sawtooth, 
whereas Cb switched between redshifts and blueshifts. Further, we noted that the Doppler 
shifts for both components oscillate during the sawtooth with an amplitude of ±20 km s“^ 
(averaging ~50 km s“^ for Cr and ~0 km s“^ for Cr) and a period of 100-200 s. This Doppler 
oscillation was found to be correlated with the spatial sawtooth oscillation, and the phase 
between these was found to be 180°. We also used the flare loops, identified in SDO/AIA 171 
A, to find the approximate conjugate point in the ER. Using an artificial slit constructed from 
the IRIS SJI 1400 A and SDO/AIA 1600 A data, we were able to identify two bands with 
the same amplitude, period, and phase as the WR sawtooth substructure. We believe the 
conjugacy of the sawtooth oscillation (along with the continuous and coherent evolution of 
the sawtooth) strongly suggests that the sawtooth is not the result of held-aligned transport 
effects, and that the two ribbons are being influenced by a common coronal source for the 
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oscillation. We use 180° phase difference between the Doppler shift and sawtooth position, 
as well as the LOS inclination of the magnetic field, to eliminate two simple wave modes as 
the source for the sawtooth. Finally, we propose an elliptical oscillation motion of the loop, 
driven by an instability in the coronal current sheet, as the mechanism behind the observed 
sawtooth ribbon substructure. 


The scenario depicte d in Figure fT^ is reminis cent of the kinds of loop oscillations some¬ 
times observed in flares (lAschwanden et al.lll999l) . Oscilla tion of loops in 171 A have pre¬ 
viously been observed in association with KH instabilities flOfman fc Thompson! 1201 ih . We 
failed to observe any oscillations in the AIA 171 A image sequence which includes Figure |3l 
but this may be due to the ~30 min delay between the loop formation, by reconnection, and 
its appearance at that cool wavelength. 


Oscillations triggered by flares are interpreted as standing MHD waves, and can have 
periods in the ~200 s range, consistent with the Alfven transit time in the long, high, weak 
loops which are often found to oscillate. This interpretation is more problematic for our case 
with short, low, strong loops. The loops, visible in 171 A (see Figure [2D appear to fit field 
lines of a const ant-alpha field with a ~ —5 x 10“®m“^, extrapolated form the line-of-sight 
HMI magnetogram. These field lines range in length from L = 40 to 60 Mm, and have field 
strengths falling to B ^ 150 G at their apices, around z ^ 10 Mm. For a typical density of 
rie — 10® cm, the coronal Alfven speed would be no smaller than 10 Mm s“^, for which the 
end-to-end transit times would be 4-6 seconds on these loops. The 140 s sawtooth evident 
at the ribbon could not, therefore, correspond to a standing wave in loops like those formed 
after the flare. 


We believe a more likely scenario to be that the short, strong field lines are moving quasi- 
statically (far slower than the Alfven speed) in response to motion imposed by an elliptical 
wave occurring at or near the current sheet. This could be the hydrodynamic motion of a 
KH instability driven by velocity shear across the current sheet fjOfman fc Thompsonll201ll) . 
For this to be the case, however, the magnetic field would have to be very nearly anti-parallel 
across the sheet, in order that a guide-field component not stabilize the instability. If the 
pattern speed, Vgt ~ 15 km/s were related to this component of the Alfven speed the angle 
between the fields would have to be within one milliradian (0.006°) of perfectly anti-parallel. 
A more likely explanation is the tearing mode (TM), which is expected to occur at current 
sheets, and also exhibits elliptical flow patterns (i.e. the velocity stream function consists of 
islands). The observed ribbon motion would be related to the inflows and outflows from the 
X-lines, estimated to be smaller than the Alfven speed by a factor, where S' 3> 1 is 

the Lundquist number of the current sheet. Flow speeds of ~ 20 km/sec would arise from 
a Lundquist number, S ~ 10®, but this may reflect the turbulent state of the sheet itself. 
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Finally, this instability would need to propagate at a speed lower than Alfven speed, but 
comparable to the flow speeds of the instability. However, we have no insights as to what 
might give the instability the precise phase velocity that we observe. 

Another point concerning the proposed scenario in Figure [12] which we left unspecified 
is the location of the loop anchor. Presumably, this anchor point is located deep enough in 
the atmosphere to allow for the observed sawtooth amplitude. We might expect that the 
Si IV, O IV, and C ll spectral lines , which are normally formed at different heights within 
the chromosphere and TR and would thus exhibit different amplitudes of the loop motion, 
providing us with some information about the anchor location. However, as we noted, these 
lines instead all display identical amplitudes of the sawtooth pattern. We speculate that this 
is due to the TR being compressed during the flare, such that the normally separated spectral 
line formation heights are all approximately the same within the ribbon, which eliminates 
any information about the location of the loop anchor. 


We also note that the small mean Doppler shift oi Cb indicate that the plasma gener¬ 
ating that component of the Si iv emission is not probably undergoing intense evaporation 
during the sawtooth. Further, the non-therma l widths for Cb are consistent with those for 
quiescent AR plasma fjPe Pontieu et al.ll2015[l . On the other hand, both components of Si 
IV show strong increases in intensity during the sawtooth, indicating the the plasma has 
been energized into enhanced emission. We do not know what is causing this brightening 
of non-evaporating plasma within the context of our proposed scenario, although we have 
speculated above that wave motion may contribute to the plasma energization. 


In their recent work, iLi fc Zhand (1201511 report an observation of ribbon substructure 
in a different flare (2014 September 10, XI.6-class) that is in some ways very similar to our 
observations of the 2014 April 18 flare. In particular, they observed a quasi-periodic slipping 
motion, which manifested as a series of bright knots which appear to move along the ribbon at 
~ 20 km s“^, close to our own estimate of Vst- T heir interpretation of th is slipping motion was 
as an apparent motion, of the kind discussed by lAulanier et al.l (1200611 and dubbed “slipping 
reconnection”. Unlike our scenario for the 2014 April 18 event, however, they speculate that 
this slipping reconnection may be the result of density variations in the reconnection region, 
which are dri ven by p-mqde osc illations above the sunspot. Also in marked contrast to our 
observations, ihi fc Zhand (1201511 report a relatively steady Si IV redshift, unmodulated by 
the sawtooth pattern. This is consistent with an interpretation a pattern projected from 
the current sheet onto the chromosphere along stationary magnetic held lines: i.e. slipping 
reconnection. Our observations, however, reveal clear motion of the plasma, through the 
Doppler shifts which are correlated with the sawtooth pattern. This rules out the hypothesis 
that the oscillations are merely projections and, at least in our flare, it rules out slipping 
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reconnection. 


It is possible that the discrepancies between these two observations arise pnrely from 
the different viewing angles of the different cases. The 2014 September 10 flare occnrred 
very close to disk center, so if there were horizontal plasma motions, tracking the elliptical 
loop oscillations, they wonld not give rise to line-of-sight Doppler shifts. Onr own case, the 
2014 April 18 flare, occurred 40 degrees from disk center, where horizontal motions would, 
and apparently do, produce line-of-sight Doppler shifts. 


The apparently contradictory observations presented in this paper and in iLi fc Zhang 


(120151) indicate that additional observations of flare ribbons, located at varying locations 
and viewing angles, will be necessary to differentiate between the two hypothesis: apparent 
motion due to slipping reconnection vs elliptical loop oscillations driven by KH or TM in¬ 
stabilities. Of course, it has only been recently that instruments with sufficient spatial and 
temporal resolution have allowed for observations of ribbon substructure. As we have noted, 
lower resolution and cadence imagers such as SDO/AIA only barely allow for recognizing 
ribbon substructure at scales less than 1", and also do not provide the Doppler information 
that allows for determining phase relationships of chromospheric phenomena. The authors 
expect that as future IRIS observations develop, particularly in conjunction with a coronal 
spectrograph such as Hinode/ElS, we will be able to distinguish between various hypothesis 
for generating ribbon substructure, and use observations of ribbons to eliminate or support 
models for reconnection in the corona. 


IRIS is a NASA small explorer mission developed and operated by LMSAL with mission 
operations executed at NASA Ames Research center and major contributions to downlink 
communications funded by the Norwegian Space Center (NSC, Norway) through an ESA 
PRODEX contract. The authors would like to thank the anonymous referee who provided 
suggestions for signihcant improvement of the original manuscript. The authors would also 
like to thank Prof. Charles Kankelborg, Dr. Sarah Jaeggli, Dr. Ying Li, Dr. Paola Testa, 
and Dr. Bart De Pontieu for helpful and productive discussions concerning the IRIS instru¬ 
ments and data analysis. This work was supported by contract 8100002702 from Lockheed 
Martin to Montana State University, a Montana Space Grant Consortium graduate student 
fellowship, and by NASA through HSR. 
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AIA 1600 2014-04-18T12:46:40.12 
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Fig. 1.— SDO/AIA 1600 A image of the flare ribbons, in reversed log^Q black-white (RLBW). 
The black inset box is the positioning of the IRIS SJI FOV, and the vertical black dashed 
line is the IRIS SG slit position. The east and west ribbons are indicated by “ER” and 
“WR”, respectively. The black “F’-shaped line is an artihcial slit discnssed in Section [331 
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HMI 2014-04-18T12:58:19.50 



Fig. 2.— SDO/HMI magnetogram of the flare active region, with standard colors. The balck 
inset box is the positioning of the IRIS SJI FOV, and the vertical black dashed line is the 
IRIS SG slit position. The black and white solid lines trace the west (WR) and east (ER) 
ribbons from Figure [T] 
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AIA 171 2014-04-18X13:12:37.47 
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Fig. 3.— SDO/AIA 171 A image of the post-flare loops and ribbons, in RLBW. The black 
inset box is the positioning of the IRIS SJI FOV, and the vertical black dashed line is the 
IRIS SG slit position. The solid white lines trace the position of the SDO/AIA 1600 ribbons 
from Fignre[Il and the black “F’-shaped line is an artihcial slit discnssed in Section 13751 
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Fig. 4.— An image of the flare ribbons from the IRIS 1400 A SJI, in RLBW, corresponding 
to the black inset boxes for Figures [T]121 The east and west ribbons are indicated by “ER” 
and “WR”, respectively. The vertical black dashed line is the IRIS SG slit, and the black 
inset box outlines where the WR crosses the slit and the positioning of the frames for Figure 
[5l The black “I”-shaped line is an artihcial slit discussed in Section [37^ 
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Fig. 5.— Time series of SJI 1400 A images, for the inset box in Figure HI showing the 
evolution of the flare ribbon in minute intervals beginning at 12:46:34 UT, in RLBW. 
Time for each image is displayed at upper left for each frame. Note that time runs from 
left-to-right in the first row, then right-to-left in the second row, and continues to alternate 
down the rows, as indicated by the arrows. 
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Si IV 1403 intensity 



2014 Apr 18 
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Fig. 6.— Upper panel; time-distance stackplot of the total Si iv 1403 A SG passband 
intensity, in RLBW. Time is given on the x-axis in UT, and the y-eocis is solar-Y in arcsec. 
Intensity scale is given to the right. Lower panel: reprint of the upper panel with a red 
outline indicating the position of the sawtooth pattern described in Section 13.11 Time is 
given on the x-axis in s (after 12:33:38 UT), and the y-axis is unchanged. The blue line 
indicates the sawtooth centroid position, and the orange line is a linear fit to the blue line. 
The numbers 1-6 indicate six peaks in the sawtooth oscillation. 
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Fig. 7.— Selected Si IV spectral lines for 28 pixels in or near the sawtooth pattern shown 
in Fignre |6l An ontline of the sawtooth is shown at npper right, and the hve rows of plots 
“A”-“E” each sample a given spatial position at several different times, as indicated by the 
“+” symbols at upper right. Each plot displays the spectrograph data for Si IV 1394 A (as 
asterisks) and 1403 A (as squares), with the time printed to the right of each plot. The 
x-axis is in km s“^ from nominal line center (same for both lines), and the i/-axis scaling is 
arbitrary (the 1394 A data has been scaled by a factor of 0.5 for direct comparison to the 
1403 A line). Additionally, the solid orange (green) lines are the two component Gaussian 
hts for the 1394 A (1403 A) spectral line, as described in Section 13.21 The dashed lines are 
the individual components Cb and Cr for each of the Si IV spectral lines, with cyan/violet 
used for the respective 1394 A spectral components and blue/red used for the respective 
1403 A spectral components. 
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Fig. 8.— Tme-distance stackplot of Doppler velocity for components Cb (upper panel) and 
Cr (lower panel) of the Si iv 1403 Aspectral line. Time is given on the a:-axis in UT for the 
upper panel and in s after 12:33:38 UT for the lower panel. The |/-axis for both panels is 
solar-Y in arcsec. Velocity scale is given to the right of each panel (scaled ±40 km s“^ for 
Cb and ±80 km s“^ for Cr). Black outline indicates the position of the sawtooth pattern, 
and gray pixels indicate discarded hts (as described in Section [3^ . 
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Fig. 9.— Plot of the Doppler velocities (in km s“^) for the two Si iv 1403 A components 
for all pixels within the sawtooth as a function of time (given in s after 12:33:38 UT). The 
black squares (asterisks) are the Doppler velocities for component Cb (C'r?), and the solid 
blue (red) lines are the mean Doppler velocity for component Cb {Cb) for all positions at a 
given time. 
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AIA 1600 




IRIS SG 1403 



Fig. 10.— Upper panel: Time-distance stack plot of SDO/AIA 1600 A intensity for a vertical 
artificial slit across the ER (shown in Figures d] and 0]), as described in Section 031 with 
labels for the sawteeth “3” and “4”. Middle panel; Time-distance stack plot of IRIS SJI 
1400 A intensity for a vertical artificial slit across the ER (shown in Figures [T] and 0]), as 
described in Section 13.51 with labels for the sawteeth “3” and “4”. Lower panel: Reprint of 
the time-distance stackplot for total 1403 A intensity from Figure [6] with sawtooth labels. 
Time in all panels is in UT. 
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Fig. 11.— Upper panel: phase portrait of the mean Doppler velocity of component Cb 
(bine line from Fignre IH]) and the relative sawtooth position (orange line snbtracted from 
the bine line from Fignre [6]). Data is indicated by the sqnares, and the connecting lines are 
to gnide the reader. Time is indicated by the color of the sqnares and connecting lines, and 
is indicated by the color bar at right. The solid black line is a linear £t to all data. Lower 
panel; same, bnt with only the data for sawtooth “4” colored (all other data in gray). The 
dash-dotted line is a linear £t to the sawtooth “4” data points. 
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Fig. 12.— Upper panel; schematic diagram of our proposed scenario, with a KH or TM 
instability in the coronal current sheet resulting in elliptical wave oscillations in the recon¬ 
nected flare loops. Lower panel: schematic diagram of the flare loop footpoint, as described 
in Section 01 showing how the elliptical wave motion relative to the LOS generates the 
observed 180° phase difference between the ribbon motion and Doppler velocities. 









